Multimode circuit QED with hybrid metamaterial transmission hnes 
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Quantum transmission lines are a central to superconducting and hybrid quantum computing[l|]. 
Parallel to these developments are those of left-handed meta-materials. They have a wide variety 
of applications in photonics from the microwave to the visible range such as invisibility cloaks and 
perfect flat lenses [2, [3| • For classical guided microwaves, left-handed transmission lines have been 
proposed and studied on the macroscopic scale. We combine these ideas in presenting a 
left-handed/right-handed hybrid transmission line for applications in quantum optics on a chip. 
The resulting system allows circuit QED to reach a new regime: multi-mode ultra-strong coupling. 
Out of the many potential applications of this novel device, we discuss two; the preparation of 
multipartite entangled states and its use as a quantum simulator for the spin-boson modelQ] where 
a quantum phase transition is reached up to finite size-effects. 
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I. INTRODUCTION 

Quantum optics addresses the interaction of quanta 
of matter — atoms — with quanta of electromagnetic 
fields — photons. This is beautifully realized in cavity 
quantum electrodynamics (QED), where the interaction 
between those units is made strong by confining the field 
into a small mode volume 0- Circuit QED takes this 
further by confining microwave photons in a quasi ID 
strip-line cavity and using supcrconductingqubits as ar- 
tificial atoms with a large dipolc moment |l|, |8| . Next to 
being a promising architecture for quantum computing, 
a multitude of basic quantum optical effects has been 
demonstrated 0] . Going beyond what can be reached in 
atomic systems, an ultrastrong coupling regime — where 
the coupling strength becomes comparable to the atomic 
energy scales — has been proposed and achieved 
[m . Furthermore in the circuit QED approach, ele- 
ments are entirely human-made and can hence be flexibly 
eng ineered. This can lead to coupling to multiple modes 
|13l - [l7t either in the same or distinct cavities. There is 
a wealth of proposals exploiting these features to create 
complex photonic states T8l-[2(i| involving a large number 
of cavities. 



II. THE SYSTEM 

In one-dimension, left-handedness is defined as the 
wave vector k; and the Poynting vector having opposite 
orientation; the phase and group-velocity are opposite. 
The triplet of electric field, magnetic field, and wave vec- 
tor, (E, H, k;), is left-handed yielding a fallingdispcrsion 
relation du){k)/dk < 0. This can be achicvcdlJ] by a dis- 
crete array of series capacitors and parallel inductors to 
ground, see Fig. [TJ A low loss left-handed transmis- 
sion line (LHTL) can be realized with superconductors 
0, mi • This is the dual (inductors and capacitors inter- 



changed) of the usual [22| discrete representation of the 
right-handed transmission line (RHTL). In practice, the 
LHTL remains a metamaterial composed of discrete el- 
ements, whereas the RHTL is a metal stri p re presented 
as the continuum limit of a ladder network [22| . Unusual 
physics arises when right- and left-handed media are in- 
terfaced We realise this with a coupled transmission 
line (CTL) shown in Fig. [U a discrete LHTL coupled to a 
RHTL to be taken into the continuum limit. This model 
is studied using Kirchhoff laws with a plane wave ansatz, 
details are in the supplementary material. Individually, 
each component has the dispersion relation 
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frequencies of the 



and Li/r are capaci- 



where bjj./i are the angular 
RHTL/LHTL respectively. Q/r 

tances and inductances as defined in Fig. [T] and Cr and 
Ir are the capacitance and inductance per unit length in 
the RHTL. Ax is the size of a unit cell. A key unusual 
feature of the LHTL, as compared to a regular RHTL, 
is the divergence of the density of modes (DoM) at a 
low-frequency bound ujm = l/2\/CiLi, seen in Fig. 2(a)[ 
implying the existence of a quasi-continuous band even 
in a cavity. In the LHTL, low frequencies correspond to 
short wavelengths due to the falling dispersion relation 
Lj{k). Thus, by only a small change in frequency, a new 
orthogonal mode can be found that is different by one 
node in the left-handed component. As the wavelength 
approaches the lattice constant, the dispersion relation 
in equation ([T]) becomes flat due to Bragg reflcction[23| 
— the DoM develops a van-Hove-type singularity setting 
the aforementioned divergence at wir. Due to the hybrid 
nature of this new CTL, the closely spaced frequencies at 
this lower band-edge have nearly-identical spatial struc- 
tures in the RHTL. The fast oscillation in the LHTL 




FIG. 1. Discrete LHTL coupled to a continuous RHTL . The regular right-handed part of the transmission line is on the right, 
connected to a left-handed line shown on the left. The terminating capacitors allow to externally access the modes. The qubit 
couples to the right-handed component. In the supplementary material it is shown that stray capacitance to ground in the left 
handed line does not change the physics. Additional stray inductances and capacitances can be taken into account as shown 
by Eleftheriades, et al. Qj. The light blue areas indicate strip-lines and ground planes. 



ensure orthogonality between modes. Figure 2(b) shows 



three consecutive low frequency modes obtained from the 
full solution. The DoM can be approximated by the sum 
of the densities in the uncoupled lines. Thus the DoM of 
the CTL is approximately 
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Ni is the number of cells in the LHTL. The agreement 
between this prediction and the numerically obtained 
modes of the full model is excellent up to small oscil- 



lations, see Figure 2(a 



The quantum behaviour of the CTL is obtained 
through canonical quantization of the circuit in Fig. [1] 
This leads to a system of uncoupled quantum harmonic 
oscillators, each described by operators a]j,a„ acting on 
modes with frequencies ujn- A qubit described in its en- 
ergy eigenbasis by Pauli matrices (7^^^ placed close to the 
CTL will couple to mode n with strength _g„ 



H/h= -^O-z + ^ gnO-x (fln 



^w„a|,a„. (2) 



If the qubit is coupled to the RHTL, ~ 9n+i for 
low frequency modes since they have similar spatial pro- 
files in the RHTL. For a flux qubit [M HI, the mode 
dependent part of the coupling strength is given by 
I?(w„)(/„(a;))/max„{(J„(a;))}. The current /„ is aver- 
aged over the spatial extend of the qubit. Figure 2(a) 



shows that the qubit can be coupled to a wide range 
of modes. This mode structure allows the qubit to si- 
multaneously couple to multiple modes when > 1 
modes fall within a frequency interval of 2gn. We refer 
to this regime as multi-mode strong -coupling. It can be 
reached with other qubits, notably transmons [2^. How- 
ever, flux qubits allow us to reach multi-mode ultrastrong 
couvlina^ldl. [ll| — gn/^n > 0.1. This regime offers many 
new possibilities for circuit QED. 



III. APPLICATIONS 

The multimode Rabi Hamiltonian, equation allows 
us to prepare multimode entangled states. Within the 
rotating wave approximation it conserves the number of 
excitations. Exciting the qubit and placing its 1 res- 
onance frequency slightly above wir, allows the qubit exci- 
tation to distribute itself over many modes, i.e., produce 
arbitrary superpositions of the form co|l; Q)+J2n Cn|0; n). 
|0; n) indicates the qubit in the ground state, a single pho- 
ton in mode n and none in the other modes. For |1;0) 
only the qubit is excited. These states are in general 
entangled as seen from their Von Neumann entropy p7|. 
Figure [3] shows the entropy, indicating multimode en- 
tanglement, of a single excitation, starting in |1;0), that 
spread out over all modes for a time t. 

The Spin-Boson model Q is a fundamental model 
of quantum dissipation which allows to understand the 
transition between coherent and incoherent behaviour as 
well as a quantum phase transition suppressing quantum 
tunnelling. It is described by the Hamiltonian in equation 
(HI in the limit where the modes form a continuum. The 
dense modes at the low-frequency end provide a generic 
and realizable quantum simulator for this model. Our 
unusual density of modes provides a novel regime of sub- 
subhomic models with a low- frequency cutoff, i.e., a spec- 
tral density of the form 
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The ground and lowest excited states are well approxi- 
mated by a multimode Schrodinger cat state of the qubit 
dressed by coherent photonic states 



|±) = i=(|L)(g)|A„)±|i?)(g)|-A„)j. (3) 



\L, R) are the eigenstates of a-^. The renormalized energy 
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FIG. 2. (a) Example of density of modes, showing a lower 
band-edge at a;iR.. Dots indicate actual modes, the blue curve 
the approximate formula which is in excellent agreement. The 
red curve shows the coupling strength between a 0.5 mm 
long flux qubit placed at a current anti-node of the trans- 
mission line, as in Fig. [1] Designing the qubit to couple 
to the 4.751 GHz, mode with strength 200 MHz, results in 
ultrastrong-coupling to 36 modes within a ±200 MHz range, 
(b) Example of the first three consecutive normal modes, the 
voltage profiles in the RHTL are almost identical. The LHTL 
has 200 unit cells measuring 100 fim each. In the LHTL the 
voltages at the discrete unit cells, for the first mode only, 
are shown by the black dots; the continuous lines serve as 
a guide to the eye to see the mode structure. Requiering 
a 50 f2 impedance and an IR cutoff at cjiR,/27r = 4 GHz 
sets C; = 398 fF and L; = 995 pH. The parameters for 
the RHTL were chosen so that it supports a full wavelength 
at ojiR. This sets the values for its total inductance and ca- 
pacitance. Therefore a 3 cm long RHTL requires a capac- 
itance and inductance per unit length of Cr — 1667 fF //im 
and Ir = 4167 pH/^m. 



splitting is 



Aeff = 
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In order for the multimodc cat state in equation ^ to 
be a good approximation to the true ground state we 
apply adiabatic renormalization Q. This suggests that 
fast modes, those with cl>„ > Acff,, adiabatically fol- 
low the qubit. Slow modes remain unaffected. Thus 

2 

An = £t0(w,i — Aeff) which leads to a self-consistency 

relation for Aoff. The ratio Aes/Ao measures the accu- 
mulated phase space distance of the dressing clouds, i.e., 
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FIG. 3. Von Neumann entropy as function of the traced out 
mode n. Here the qubit is put in the bath at a frequency ujq > 
cJiR and the system left to evolve for a dimensionless time tg. 
The Von Neumann entropy is computed for the system after 
having traced out the qubit and mode n. We then vary n 
to see how much the entropy increases from Eq, the entropy 
when only the qubit is traced out. Eq > indicates that 
their is at least bi-partite entanglement. When additionally 
tracing out mode n the entropy E„ increase above Eg and 
this for all modes. The latter indicates complex multi-partite 
entanglement. 



the total cat size Q by taking the logarithm of equa- 
tion 11]). Thus, the low-energy states of the system are 
strongly renormalized as are the effective energy scales. 
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FIG. 4. Plot of the renormalized tunnelling as function of the 
global coupling strength g. The individual mode coupling g„ 
is the product of g and the mode-dependent spatial contri- 
bution bounded by unity. This allows simultaneous tuning 
of the coupling to all modes. The discontinuous drop in Actt 
distinguishes the regions of weakly and strongly renormalized 
energy splitting. 

A true dissipative quantum phase transition 0, [l^ has 
Aeff = in the localised phase. This limit would be 
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reached if the modes were infinitely close (hence arbitrar- 
ily close to wir) as would result from an infinitely long 
LHTL or if wir — )■ as in the case of infinitely dense 
LHTL unit cells. Note that in the usual sub-Ohmic spin- 
boson model, the latter is assumed. We thus conclude 
that our system approaches a quantum phase transition 
in the infinite sample limit. 

To corroborate the finite size-behaviour, we have stud- 
ied the ground and first excited state of the qubit-CTL 
model using its actual modes in the adiabatic renormal- 
ization approach. We identify multiple regimes: for weak 
coupling or large Aq/wir, there is only weak dressing 
manifest by a small shift of Aq. At stronger coupling, 
we observe the quasi- localized phase, with Aeff <^ Aq. 
Remarkably, even at finite length, the two regimes are 
separated by a discontinuous transition as indicated by 
Figure SI Figure [S] shows the corresponding finite-size 
phase diagram highlighting the need for ultrastrong cou- 
pling. We see that by tuning the bare qubit frequency 
slightly above the cutoff, Aq > wir [1^ we can tune the 
system through the phase transition in situ. 
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IV. CONCLUSION 

Wc have proposed an engineered hybrid transmission 
line that allows to reach a new multimode strong coupling 
regime of circuit QED by combining a regular line with 
a metamaterial. This will open the way for novel appli- 
cations in microwave photonics and strongly correlated 
photon states, out of which we have outlined the genera- 
tion of multimode entanglement, multimode Schrodinger 
cat states and quantum phase transitions. 



V. ACKNOWLEDGEMENTS 

We acknowledge useful discussions with Heiko Rieger 
and Britton Plourde. We thank Emily Pritchett for her 
careful reading of the manuscript. This work was funded 
in parts by DARPA through the QuEST and lARPA 
through the MQCO program. 



Appendix A: Electric Properties of Right and 
Left-Handed Transmission lines 



The following section presents the electric properties of 
the composite transmission line. The derivations follow 
typical microwave engineering as outlined in Chapter 2 of 
Ref. m. We show how the system behaves at the bound- 
ary, i.e. under what conditions energy flow across the 
boundary is optimized. Examining the relation between 
voltage and current in each medium gives the impedance 
of the lines; the latter are then used to express the reflec- 
tion and transmission coefficients at the boundary. They 
are found by imposing continuity of voltage and current. 
Then the voltage profile along the CTL is discussed. The 
section finishes in showing that their is a dense mode 
spectrum at a low-frequency bound. 



FIG. 5. Finite-size phase diagram of the qubit-CTL system. 
The density of states for the case at hand is the same is in 
Fig.[2R| 



On the level of partition functions, this model is equiv- 
alent to a one-dimensional Ising chain [2^, [sO] . This is 
discussed in more detail in the supplementary material. 
In the present case, this would be an Ising chain with 
an interaction that decays as cx |z — jl^^^^, where i and 



j are site indices, up to a range r cx , after which 
it decays exponentially. Thus, when cooled from high- 
temperatures the system is well described by mean-field 
theory, which predicts a ferromagnetic phase transition, 
until the correlation length reaches r. At that point, the 
system follows short-range physics and remains paramag- 
netic between magnetized blocks of size r — in analogy 
to the tunnel coupling in the spin-boson model falling 
deeply, but not to zero. 



1. Equations of Motion 

The right and left-handed parts of the transmission 
line are respectively shown in the right and left side of 
Figure [1] of the main text. In terms of magnetic flux $ 
their uncoupled Lagrangians [H, HH are 



n<0 



ri>0 



Cr Lr are the capacitance and inductance of the dis- 
crete model of the RHTL which will then be taken to the 
continuum limit. Ci and Li are the capacitance and in- 
ductance of a unit cell of the LHTL. Through the Eulcr- 
Lagrange equation, these Lagrangians produce the fol- 
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lowing equations of motion 
1 



ClLl ( $„_1 - 2$„ + $„+l ) = 



(Al) 
(A2) 



The latter are none other than Kirchoff 's current law for 
each cell of the transmission line. This is seen by noting 
that, for the case of the RHTL, the voltage on capacitor 
n is Vn = — <&n and the current flowing through this 
capacitor is /c„ = CrVn- The flux difference between two 
nodes relates to current through $„ — = Lrln-i- 

Thus equation (|A1[) is none other than /„ — = /c„i 
see also figs. 6(a)|6(b)" Similar reasoning applies to the 
LHTL. 



In-l 



n - 1 n 



-fc„_i = CrVn 



(a)RHTL unit cell 





(b)LHTL unit cell 



FIG. 6. Unit cells of of both transmission line with Kirchoff 's 
current law. 



2. Dispersion Relation and Wavelength 



3. Brillouin Zone of the Transmission Lines 

Given that the transmission line is made of discrete el- 
ements of length Ax the minimum distinguishable wave- 
length is Amin = 2/S.x. In the case where both lines are 
made of discrete elements their is both an IR and UV 
cutoff. In the case where the RHTL is continuous, the 
UV cutoff is set by the superconducting gap of the ma- 
terial, this shall be ignored given that it happens at fre- 
quency scales that are not relevant when coupling qubits 
to transmission lines. Measuring distance in units of the 
unit cells, i.e. Ax = 1, gives the maximum allowed wave 
number as 
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In the left-handed medium an upper bound on the max- 
imum wave number corresponds to a lower bound on the 
frequency since fc; ^ f/w/. In the right-handed medium 
A^max gives an upper bound on the frequency 
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Therefore due to the discreteness of the composite TL, 
the allowed modes fall in the interval [wirjOJuv] — the 
Brillouin zone of the discrete CTL. In what follows, it 
shall be assumed that the RHTL is continuous and has no 
upper frequency bound other than the superconducting 
gap. 



The dispersion relations in the discrete RHTL and 
LHTL are found from the equations of motion [XT] bv as- 
suming a plane wave ansatz $n(i) = exp{— i(fcnAx— wi)} 



and using the relation 2— e 
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Ax is the length of a unit cell. uJr/i is the angu- 
lar frequency of a wave with wavcvector k^/i in the 
RHTL/LHTL. In the right-handed medium, the usual 
continuous dispersion relations vk ~ lo can be recovered 
by taking the limit Ax 0. With the inductance and 
capacitance per unit-length Z,. = L^/Ax and Cr = Cr/^x 
respectively, the square of equation (jA3p becomes 



4. Current, Voltage and Impedance 

In the following section the relation between cur- 
rent and voltage is examined in order to find the 
impedance. Making use of Kirchoff's laws and consid- 
ering the RHTL's unit cell shown in Fig. 6(a) gives 



(A7) 



In the LH unit cell, shown in Fig. |6(b)[ current and 
voltage are related through 



1 f i 
— / VrAt = —Vn 



a. Plane Wave Ansatz 



(A8) 



a;^ZrCr(Ax)^ = 2 - 2cos(fcAx) w (fcAx)^ 

Cancelling (Ax)^ on both sides and identifying v = 
Xj \JlrCr as the phase velocity of the waves gives the usual 
continuous dispersion relation for the RHTL. 



The voltage along the CTL can be decomposed into 
two parts, each originating in one side of the CTL and 
propagating towards the boundary. The waves propagat- 
ing along the CTL must have the same energy therefore 
we set = a;^ =0;. originates in the RHTL whilst 
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V comes from the LHTL. Explicitly, these two dimen- 
sionless voltages with unity as maximum amplitude are 
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When the waves encounter the interface at = they are 
partially transmitted and reflected; and are the 
reflection and transmission coefficients for V^] they will 
be determined in section iA 51 It is important to note that 
in the left-handed line the group velocity is negative: the 
sign of the wave vector is opposite to the Poynting vector. 
The voltages given by (|A9I) and (|A10I) are associated with 
currents that take the form 
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b. Transmission Line Impedance 

Following mi the impedances of the RHTL and LHTL, 
Zr and Zi respectively, are 
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When taking the limit, the dispersion relation was used 
and the usual form for the impedance of a continuous 
strip line is recovered. Here the definition of wir, i.e. 
equation (|A5[) was used. This shows that the quantity 
under the square root is always positive. Furthermore 
the impedances enjoy the property that = y^Lr/Cr 
and \Zi\ = ^ Li/Ci. The terms featuring a frequency de- 
pendency are the result of the discreteness of the model. 
This can further be confirmed by taking the long wave- 
length limit; when A ^ Ax it is expected that the dis- 
crete transmission lines behave like continuous ones. For 
the LHTL the long waveleng ths are achieved by the limit 
w oo for which Zi — >■ yjLi/Ci. Making use of the 
impedance, the currents in relations (jAlip and (|A12p 
can be more conveniently written as 



The six coefficients and Ip are introduced to re- 

late voltage to current. They are found by making use of 
Kirchoff 's laws (jA7p and (|A8p . For instance as an explicit 
example is found by considering the current difference 
between cell n and n + 1 
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Where the first equality in (jA13[) follows from the plan 
wave ansatz and the second one is a result of Kirchoff 's 
law. The six coefficients are then 
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Relations ()A14p and ()A16P were found with Kirchoff's 
current law for the RHTL and (|A15|) and (|A17|) with the 
LHTL current law, i.e. equations (|A7p and (|A8|) respec- 
tively. 



5. Reflection and Transmission Coefficients 

The exact form of the reflection and transmission coef- 
ficients is found by imposing voltage and current continu- 
ity across the boundary between the RHTL and LHTL. 
This results in 
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Solving these two systems yields the reflection coefBcients 

Zi — Z* Zr 



r+ = 



r- = 



Zi + z,. z* 



Zr — ZI Zi 

Zi + Zr ZT 



They enjoy the property |r+p = |r~p. Given that the 
CTL is a passive system energy flowing in the RHTL 
must be equal to energy flowing in the LHTL. Mak- 
ing sure of this serves as a consistency check. The 
time averaged power along a transmission line is = 
Re[y^Z^*]/2. It is important to note that we expect the 
outcome to be independent of the position along the CTL 
since their are no power sinks or sources. Thus, in terms 
of time averaged power, the following must hold 



ri<OJ 



>0^n>o] 



If the + solution is considered, i.e. voltage given by (|A9[) 
and current by (jA19|l . then the energy flows in the RHTL 
and LHTL are respectively 
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2\Zr?\Zr^Zi? 



P: 



-Re 
2 



_ 2 \Zr\^ Re {Zi) + Re + Zf)) 
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The two terms with r+e~*'^' /Zr arise from the interfer- 
ence between incoming wave and reflected wave, however 
they do not contribute to the power since A — A* Ci i M. 
By comparing these two equations it is immediately seen 
that = P^+, i.e. energy is conserved. The same rea- 
soning applies if the — solutions of voltage and current 
are considered. 



The flow of power across the boundary is thus maxi- 
mum when the impedances are matched in magnitude 
i.e. Lr/Cr = Li/Ci. In this case and since the voltage 
is taken to have unit amplitude, the maximum power is 
(2Zr)~^. This corresponds to the limit w — oo for which 
the two lines are perfectly matched. Additionally and re- 
gardless of impedance, the following properties hold when 
going closer in frequency to the IR cutoff 

lim P,,v = 



lim 



dP^ 

duj 



The power flow as function of frequency increases mono- 
tonically. The infinit slope of Pav(w) at wir guarentees 
that the first mode will have non-zero power flowing 
across the boundary. An example of the power across the 
boundary of a 50 SI CTL is shown in Fig. [71 The black 
dots represent the normal modes, discussed in section 
lA 71 The lowesst mode obtained in a real geometry, with 
frequency loi/2tt — 4.149 GHz has a power flow of 42.1% 
of the maximum despite the fact that Pa^vi^m) = 0. 
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FIG. 7. Power flow across the boundary as function of mode 
frequency. The black dots are the numerically found modes 
discussed in Fig. 2(a) of the main text. 



a. Impedance Matching 

For the two parts of the line to be able to communicate, 
the energy flow across the boundary must be greater than 
zero. Using Zf = Lr/Gr and writing P^ as function of 
frequency yields 




6. Voltage on the Composite Transmission Line 

In the following boundary conditions are imposed on 
the voltage allowing a normal mode decomposition to be 
found. The normal mode indexed by k has a frequency oJk 
and wave vectors kr^i{uik)- The voltage along the line, for 
mode k, at node n and time t is given by a superposition 
of V+ and V~ 

ik) 

The coefficients Vq and lS.(j)k are both real and are in- 
troduced so as to match the phases and amplitudes of 
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both waves. They depend on the frequency of the mode 
considered. The discussions in the last sections were re- 
stricted to the interface and in the lines but did not in- 
clude effects linked to the finite length of the CTL. Here 
shall be considered the case where the line is asymmetric; 
the node in the right-handed line furthest away from the 
interface is labeled N,. and the one for the left-handed line 



is —Ni 7^ —Nr- The boundary conditions will determine 
what modes the CTL can support. Similarly to appendix 
A in the boundary conditions for an open-circuited 
CTL are given by charge neutrality and impose that the 
spatial integral of the voltage vanish at the extremities of 
the TL. The integrals 6^ and 6*+ are respectively found 
from equations (|X9)) and (jAlOp 



iki 



ikj. 



n < 
n > 



:{t) 



n<0 



iki iki — 



At the extremity of the RHTL and LHTL 9 for mode k is respectively given by 



7N,.,k — 



-Ni,k — 



1 
1 

iki 



-iujkt 



The boundary conditions are |0-Ar^,fcp = \0Ni.k\'^ ~ 
0; implying that a solution for the following non- linear 
system of equations has to be found 



-2ikrNrAx 



Y+ _ yWrp-^iArjik 



yif^) gi2k,NiAx _ Y^'''>Y~ - J^+g-''A0fc 



(A21) 



Where it is understood that F*, T*, fc,,, fc;, V^'''"^ and 
A(f)k all depend on the mode frequency cok- This system 
has four parameters Cr.i and i^,; and three unknowns. It 
is solved numerically by selecting the desired line capac- 
itance and inductance and searching for (V^^*'' , A0fe, w^) 
triplets. As example, three modes supported by a 50 51 
CTL with Cr = 1667 iF/fim, Ir = 4167 pH ^m Ci 
398 fF and Li ^ 995 pH are shown in Fig. |2(b)| of the 
main text. The LHTL is made of 200 cells and the RHTL 
is 3 cm long. The values for the capacitance and induc- 
tance were chosen so that the IR cutoff is at 4 GHz. The 
voltage at the unit cells is indicated by black dots and 
the coloured lines are to help visualizing the modes. As 
can be seen from the figure, the composite line can sup- 
port modes with slightly different frequency. For these 
neighbouring modes the voltage maximums in the RHTL 
are only slightly shifted in position. This is the result of 
the wavelength in the LHTL being very small. A small 
change in frequency is thus enough to include and extra 
anti-node in the LHTL. For the high frequency modes 
this trend is reversed the wavelength is short in the RHTL 
and long in the LHTL. The similar voltage profile in the 



r 



RHTL at low frequencies is what will allow cQED to 
reach the multi-mode ultrastrong coupling regime. 



7. Density of Modes 

The density of modes (DOM) is found by numerically 
solving equations (jA21[) . An example of DOM for a CTL 
is shown in Fig. 



2(a) of the main text. A theoretical 



understanding can be gained by considering the right- 
handed and left-handed transmission lines separately. 
The DOM gives the number of modes within the interval 
fc to fc + dfc by I?(fc)dfc ~ dn. The frequency dependent 
DOM is 



dn 1 
dfc doj/dfc 



In the case of a chain of N identical cells of length 
L = NAx the mode numbered by n has a wave num- 
ber fc = 7i27r/i thus dn/dk = L/2tt. Using the disper- 
sion relations given by (|A3[) and (jA4p the DOM for the 
individual systems are 



Vr (w) = 
Vi (lu) = 



tan (pi sin cpi , 



ki{uj)/S.x 



Ni is the number of unit cells making up the LHTL. The 
DOM is then approximated by I?(w) = 'Dr{uj) + Viiuj). 
As shown by Fig. 2(a) of the main text there is a good 
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agreement between the approximate and the immerically 
exact DOM up to some small oseillations. Using the 
dispersion relations to make explicit the dependence on 
w yields 



2?M 



2Ni 



^IR 



dr is the length of the RHTL taken to the contiimum. In 
practice and close to the IR cutoff the last term, com- 
ing from the RHTL, can be neglected since the DoM is 
dominated by the LHTL. This is what we call a sub-sub- 
ohmie density of states. Indeed, expanding the functions 
in the DOM as a power law around cjir. Expansion in 
5uj = uj — wiR+ reveals that 'Di{5u)) cx (Jo;)^^''^. Drop- 
ping second order terms in 5uj results in 



V (wiR + 5ijj) 



Ni 



1 



account through the finite size factor defined by 



dxVu{x) 



max / AxVk {x) 

The integral of the voltage is taken over the qubit's spa- 
tial extent and normalized to the mode that couples 
strongest to the qubit. If the qubit should be placed 
at a current anti-node, such as flux qubits, then volt- 
age is replaced by current in the definition of the size 
factor. Fk shows which modes couple the most to the 
qubit. Given the mode structure a qubit coupled in the 
RHTL (LHTL) will couple to many low (high) frequency 
modes. The case of a 0.5 mm long fiux qubit coupled to 
the RHTL is plotted in Fig. 2(a) of the main text by the 



red line showing FkD {u>k)- Using this DOM and design- 
ing the flux qubit to couple to the 4.751 MHz mode with 
strength 200 MHz. results in ultrastrong coupling to 36 
neighbouring modes within a ±200 MHz range. 



Appendix B: Coupling a Qubit to the Composite 
Transmission Line 

As shown in the section lA 6| neighbouring modes 
with low frequencies have close spatial structures in the 
RHTL. Indeed in the LHTL a low frequency mode has a 
short wavelength; a small change in frequency is sufficient 
to produce an extra node. At high frequencies the situ- 
ation is reversed; the wavelength is short in the RHTL 
and long in the LHTL. Given that for neighbouring low 
frequency modes in the RHTL, the anti-nodes are close 
together, it is interesting to put a qubit close to the CTL. 
Thus coupling it to many modes. The situation of multi- 
mode ultrastrong coupling arises when their are > 1 
modes within an interval of 2gk and Ofc /j-^k > 0.1. This 
situation is realized with fiux qubits jl2| . The last con- 
dition corresponds to coupling strengths for which the 
rotating wave approximation starts breaking down. In 
some situations it can also be desired to have only strong 
coupling between the qubit in the CTL. In this ease the 
qubit will still interact with many modes but the RWA 
will hold. 



1. Qubit Size Factor 



2. Multimode Ultrastrong Coupling and the 
Spin-Boson Model 

It has been experimentally demonstrated that flux 
qubits can achieve the regime of ultrastrong coupling 
where the RWA breaks down [ll|, [l3| . A theoretical sur- 
vey of how this is achieved is described by Bourassa et 
al. . The flux qubit has a double well potential; when 
coupled to a CTL the system is described by a Spin- 
Boson model 

F[sb = 

^ ^ M M 

- -Aqt^ + -ef, + uJkola^ ^^k + h) 

k=l k=l 

CLk and are the creation and annihilation operators for 
mode k of the CTL. Here we have, compared to the con- 
vention in the rest of this work, rotated the basis for 
the two-state system according to ffz/x = t'x/z in or- 
der to connect to the spin-boson literature The coupling 
strength is written as gFk to be able to simultane- 
ously tune the coupling to all the modes whilst keeping 
the effect of the qubit's size. In the following the sym- 
metric case shall be assumed by setting the bias to zero, 
i.e e = 0. When coupling the qubit to the modes close to 
the lower frequency limit ojir the density of states is well 
below the sub-ohmic limit (in this case s ^ —1/2 < 0). 
The coupling with the bosonic bath leads to a renormal- 
ization of the tunnelling element Aq i— > Aeff. 



To couple a qubit to as many modes as possible, 
the anti-nodes must move as little as possible. This is 
achieved in the RHTL at low frequencies close to wir. 
Thus if Ag is the qubit's transition frequency, loyr, ^ Aq 
is desired. The effect of the qubit's size can be taken into 



a. Approaching a quantum phase transition 

Showing that the qubit-CTL system approaches a 
quantum phase transition up to finite size effects, is done 
by making use of the map between the Spin-Boson model 
and the ID Ising chain [30|. A true quantum phase tran- 
sition i.e. Aeff = is reached only in an infinitely dense 
LHTL. Following the reasoning in Weiss [2^ these parti- 
tion functions are expressed as sums over the number m 
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^^l^=EE (l^j /^-'^m / d,S2„-l... / dsi XCXp<^^^(-ir+^>(s,-5,) 
±90 m=0 



2rn j-1 







{j = 2 i=l 

2m j-1 



pi PA PA I J ^ 

m=0 n n n I J=2 1=1 



In the flux qubit-CTL system a kink corresponds to 
the system tunnelling from one side of the potential well 
to the other. In the ID Ising chain a kink corresponds 
to a spin domain wall. By comparing the two partition 
functions, it is seen that the second integral of the Ising 
interaction potential U plays that same role as the kink- 
kink interaction of the Spin-Boson model W{r). There- 
fore the quantity d^W{T) is equivalent to the interaction 
potential of the Ising model. Using the definition of W 
gives 

oo 

V{t) = dlW{T) = adl J dio'^ (1 - e--^) 



aN, 



1 - e- 



V2wiR7r J uj'^^uj - ujiji 



aN] pK 



\/2t^iR7r 

The constant a includes all the dependencies of qubit- 
CTL coupling. An example of which can be found in 
[lo| . Here the density of states was approximated by the 
w"^/^ law valid close to the IR cutoff. The interaction 
potential has an exponential cutoff, and it is know that 
this model has no phase transition. In the Ising model 
this corresponds to forming domains of size 1/wir. These 
domains then only posses nearest neighbour interaction, 
for which their is no phase transition. In the case of a flux 
qubit coupled to the CTL, the bare tunnelling element 
Ao is renormalized until it is pushed below the IR cutoff. 
At this point, their are no more modes left to further 
reduce Aoff. 



6. Adiabatic Renormalization 

It is expected that if the qubit is tuned so that 
Ao > LJm the energy splitting will get renormalized to 
Aoff. This process can be well approximated by adiabatic 
renormalization This process abruptly stops when 
Aeff < wiR. It is expected that the system is described 
by the state 



l±) 



1 



V2 



^|A.)±|1)(^|-Afe) 
k fc y 



I 

where the coherent states \Xk) and |— Afc) are created by 
the displacement operator I?(Afc) = cxpjA/ca^ -f- A^dfc} 
with a sign depending on whether the qubit is in the 
|0) or |1) state. Adiabatic renormalization suggests that 
Xk — glQ{uj — Acfi)/w^.. The set of modes with uj^ > Ao 
renormalize the tunnelling element 



A' = Aoexp I -2 
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This results in a new set of modes satisfying Aq > Wfc > 
A'. The latter formula can be iterated to find the fully 
renormalized tunnelling element Aeff. Figure HI shows the 
renormalized tunnelling element as function of the cou- 
pling strength. The finite size of the qubit was also taken 
into account according to section IB II As can be seen 
there is a sudden drop in Aoff once some critical coupling 
strength ^crit is reached. The renormalized coupling cl- 
ement does not drop to zero but to some finite value as 
expected from Section IB 2 al The values of Ocrit corre- 
spond to what is experimentally achievable jl2l | indeed, 
Nicmczyk et al. have measured g/l-K = 636 MHz with 
a;/27r = 5.357 GHz. By identifying the regions where 
Aoff < wiR in the (Ao,(7) plane shows regions where this 
strong renormalization occurs. The border between the 
two regions is accompanied by a discrete jump in Aoff as 
shown by Fig. |4]of the main text. 



Appendix C: Physical Realization 

For the CTL to reach the quantum regime it should 
be superconducting [H . Here we propose to fabricate the 
coupled transmission lines by standard fabrication of the 
RHTL and using Josephson Junction Arrays (JJA) to 
create the LHTL as Josephson junctions provide a com- 
pact way to engineer sizeable inductances [32| - [33 | . These 
Josephson junctions are large enough to be used in the 
linear regime only. The values of the critical current and 
capacitance are determined by impedance matching be- 
tween the two parts of the line to ensure maximum trans- 
mission, and the desired IR cutoff. The impedance of the 
RHTL is Zr = \Jlrlcr and should be chosen to be as 
close as possible to Zo = 50 51 to match the impedance 
of the electronics. For generating multi-mode entangled 
states and studying the Spin-Boson model, the IR cut- 
off frequency should be close to the qubit transition fre- 
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qucncy Ag and the impedance close to Zq. Thus 
using impedance matching and the definition of wir, i.e. 
equation (jA5[) . the parameters of the LHTL unit ceh are 



Ci 



1 



and L, 



2a;iR 



Most superconducting qubit transition frequencies are in 
the range 4 to 12 GHz placing the capacitance in the 
range C; € 398 fF to 133 fF. The inductance should be 
from 995 pH to 332 pH respectively. To first order the 
inductance of a JJ is 



2e/o 



A desired inductance in the range 995 pH to 332 pH 
corresponds roughly to a critical current in the range Ic G 
330 to 990 nA. For a fixed cross section and material, 
the length of the RHTL has the effect of changing the 
fundamental frequency. For the continuous RHTL alone 
the frequencies of the normal modes are Uk — k-K j \fcrU-dr 
where and are the capacitance and inductance per 
unit length respectively. As an example it can be desired 
to have a voltage anti-node in the middle of the RHTL. 
This sets the capacitance and inductance per unit length 
to 



27r 



LOlB.Zodr 



and Ir 



2ttZo 

LUmdr 



An IR cutoff at 4 GHz and a 3 cm long RHTL sets c,. = 
1667 fF//im and Ir = 4167 pH//im. 



The Josephson Junctions making up the inductors of 
the LHTL additionally have a capacitance in parallel to 
the inductor. The unit cell is thus more properly mod- 
elled by the one shown in Fig. [51 Taking this capacitor 
into account, the Lagrangian of the LHTL now reads 



n<0 



The equation of motion for which is 



The ansatz $„(t) ~ $o exp{— i(fc;nAa; — ujt)} produces 
the dispersion relation 



1. Parasitic Capacitances due to the JJ's 



n-1 




FIG. 8. Modified left-handed unit cell taking into account the 
capacitance of the JJ. 



1 



kiAx 



4Ci 



-1/2 



At the infra-red cutoff sin(fcA/2) = 1 thus as long as 
Cjj ^ iCi the infiuence of the JJ's capacitor will be 
negligible. 
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